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The duality relation for the effective conductivity cr e of 2D isotropic hcterophase 
systems is used for obtaining the exact results for cr e at arbitrary number of 
phases N. The exact values of <r e correspond to the fixed points of the duality 
transformations. The new exact results for <r e , generalizing the well-known 
exact values of cr e for N = 2, 3, are found. It is shown that for N > 3 there 
exist the whole hyperplanes in the space of phase concentrations, on which <r e 
takes constant values. These results are checked in the framework of various 
approximations for different random heterophase systems. 
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1. Introduction 

The properties of the electrical transport of classical inhomogeneous (ran- 
domly or regularly) heterophase systems, consisting of N (N > 2) phases with 
different conductivities at (i = 1,2, ...,N), have a great significance for prac- 
tice. The main problem in theory of the electrical transport in these systems 
is a calculation of the effective conductivity a e at arbitrary conductivities and 
phase concentrations yQ. Despite of seeming simplicity of this problem only a 
few exact results have been obtained so far and all of them for two-dimensional 
systems pOE]- They are based on the exact duality relation, which is a conse- 
quence of the well-known duality between potential and tangent (or divergent 
free) fields in two dimensions. This very important relation has been obtained 
firstly for isotropic two-phase systems, where it can be represented in the form, 
connecting cr e at adjoint concentrations x and 1 — x [21 EI 

a e {x,ai,a 2 )cr e (l - x,ai,a 2 ) = cria 2 . (1) 

It allows to find the exact universal, remarkably simple, formula for cr e at equal 
phase concentrations X\ = x% = 1/2 

c e = \fo\02. (2) 

Analogous exact result for er e , coinciding with (2), has been obtained in N = 3 
case by Dykhne when o\02 = 03 and the concentrations of two first phases 
are equal x\ — xi — x. These exact values of <r e turned out to be very useful 
for various applications of the theory, especially the equation (2), which has 
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been implemented for checking of various approximate expressions for <j e and 
in the percolation theory The dual relation (1) admits the generalizations 

on anisotropic systems and polycrystals, systems in magnetic field El E]- but 
it cannot give exact results for them (except isotropic systems in magnetic field 

W 

In this letter, using a general duality relation, we will obtain exact results 
for <r e of 2D isotropic heterophase systems with arbitrary number N of phases. 
They generalize the above results for N — 2,3. In particular, it will be shown 
that for N > 2, contrary to the N — 2 case, there exist the whole hyperplanes 
in the space of concentrations, on which a e takes constant values. The obtained 
results are checked in the framework of various approximate schemes. 

2. Symmetry, duality and fixed points 

We begin with a discussion of some general properties of the effective con- 
ductivity <r e of N -phase systems with the partial conductivities Ui > (i — 
1,2, ...,N) and concentrations Xi, which satisfy the normalization condition 
SiLi x i — 1- The effective conductivity must be a symmetric function of pairs 
of arguments ( Xi, o\ ) and a homogeneous (a degree 1) function of <7j. For this 
reason it is invariant under all permutations 

a e (xi,ai\x 2 ,cT2\...\x N ,a N ) = a e (P io (x ll cr 1 \x 2 ,a 2 \...\x N ,cr N )), 

where Py is the permutation of the i-th and j-th pairs of arguments, and can 
be represented in the next form 

a e (xi,ax\x 2 ,(T2\...\x N ,a N ) = a s f(x 1 , z x \x 2 , Z 2 \.-\x n , z N ). (3) 

Here Zi = ai/a s , where <j s is some normalizing conductivity. For our purposes 
it will be convenient to choose it in the symmetrical form a s = (sn) 1 ^, where 
sn = Y\i=i a i i s the N-th elementary symmetrical function. Then 

z i = a i /{s N ) 1 ' N , if zi = l. 

As a result the function / is a symmetrical function of N pairs (Xi, z,), only 
N — 1 of them are independent. 

The effective conductivity of N -phase systems must satisfy the following 
natural boundary conditions: 1) a e of N -phase system with n (2 < n < N) 
equal partial conductivities must reduce to er e of system with N — n + 1 phases 
and the concentrations of the phases with equal conductivities must add; 2) it 
has not depend on partial <Ti and must reduce to the effective conductivity of 
(N — 1) -phase system, if the concentration of this phase x< = ; 3) it must 
reduce to some partial erj , if x% = 1. 

A duality relation for a e exists for any two-dimensional N -phase system 
(see, for example, |3IH]). Its most general form for a e valid at arbitrary phase 
concentrations is 

a e (x 1 ,a 1 \...\x N ,a N )a e (x 1 ,a 2 /a 1 \...\x N ,a 2 /(T N ) =<r 2 . (4) 
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where <jq is some arbitrary constant conductivity, characterizing a duality 
transformation. 

It is useful to note that the relation (4) can be written also in the other form 

er e (xi,(7o 2 /(Ji\--\x N ,cr 2 /a N ) = a 2 /a e (x 1 , ai\...\x N , a N ), (4') 

which means that a e at the inverse values of the partial conductivities Oi is 
equal to the inverse a e . For this reason it can be named also the inversion 
relation for effective conductivity. 

The mutually inverse values are related: their product is equal a 2 (Fig. la). 
We will use also, for brevity, the notations of the form {x} for all whole sets of 
variables like Xi (i = 1, N). 

An existence of the exact duality (or inversion) relation is very important 
property of two-dimensional inhomogencous systems. It can be used for the 
obtaining of exact values of effective conductivity. For this one needs to know 
the properties of the partial inversion transformations 

jW-.ai^al/au (J«) 2 = / (i = 1, 2..., N), 

where I is the identity operator. Each partial inversion transformation acts 
only on the partial conductivity of the corresponding phase. Strictly speaking, 
the inversion transformation is well defined only for all Ui ^ 0, but one can 
extrapolate it formally on the case m = 0. In this case it transforms a system 
with <7j = (a dielectric phase) into a system with Ui = oo (an ideal conducting 
phase). 

(i) 

The transformation J& has only one fixed point Oi = <jq. The i-th and 
j-th partial inversions, having the same parameter a 2 — OiOj (we denote them 
jty and respectively), are very important, because they change Ui into 

<jj and vice versa 

Jifoi = cri<Jj/<Ji = <Jj, j\f<jj=<Jicrj/crj=<Ji. (5) 

The physical sense has only the product of all partial inversions with the same 
parameter cr 

N 

•ko-n^o. j i= i - ( 6 ) 

Giving to <7q different values, one can obtain various relations, connecting a e 
at different values of its arguments, which can serve for checking of solutions 
found by other methods. 

Since in our problem J aa acts on a space of symmetrical homogeneous 
functions it is useful to consider its action on the basic elementary symmetrical 
functions s k {{cr}) ({a} = (a lt ...,a N )) 

N i k 

s *= E 11^ (k=l,...,N). (7) 

il<i2<---<ik h 
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It is easy to see that 

Ja Sk = (?l k S N -k/s N - (8) 

It follows from (8) that only sn transforms into itself. Since a e is a homoge- 
neous functions of degree 1, it will be convenient to pass (for oi ^ ) to the 
projective variables Zi and to introduce projective basic symmetrical functions 

Sk 

h = s k ({z}) = s k /(s N ) k / N , k = l,..,N-l, s N = l. (9) 
Then we have from (9) 

Ja Zi = 1/Zi, Ja S k = S N -k, (10) 

i.e. the inversion J ao for any cr inverses Zi and interchanges adjoint §k- This 
property strongly simplify a search of possible fixed points needed for obtaining 
the exact equations for a e . In particular, the duality relation for functions / 
has the next form 

f({x},{z})f({x},{l/z}) = l, (4') 

We use the general duality relations (4,4') for the obtaining an exact closed 
equation for a e at some special values of <7j and Xi. It follows from (4) that 
such equation can be obtained if one finds a transformation T( J, P) , consisting 
from inversions and permutations, and a combination of arguments ({#'}, {<r '}) , 
which is a fixed point of T 

T(J, P)<j e ({x'}, {a 1 }) = a e ({x'}, {a 1 }). (11) 

Then, implementing this transformation to <r e ({x'}, {&'}) in (4), one obtains 
the following exact equation for a e at these arguments 

<T 2 e ({x'},{a'})=al (12) 

which gives an exact value of a e 

<7 e = cr - (13) 

Here cto is a parameter of the inversion, entering in T(J,P). The corresponding 
exact equation for function / has the form 

f\{x'},{z'}) = l, (14) 

i.e. the exact value of / at any fixed point is equal to 1. Consequently, the 
possible exact values for a e at these fixed points must satisfy the next equalities 

a e {{x'},{a'}) = a a = {s N ) 1/N . (15) 
If T( J, P) contains the interchanging inversion then 

er e = (T = = ( s n) 1/N , i ^ 3- (15') 
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Fig.l. (a) A schematic plot of possible values of a%, (i = 1, 2), satisfying the relation 
<j\02 = erg with some fixed value <tq ; (b) a concentration phase diagram for N = 3 , 
the bisectrisses correspond to the pairwise equal phase concentrations. 



Thus, (15,15') give a general form of the exact values of a e in these fixed points. 
The last equalities in (15) and (15') give a strong constraint on possible fixed 
point values of Cj. 

In case of N — 2 such fixed point exists at arbitrary <ii (i = 1,2), but 
only at the equal phase concentrations x\ = x% = 1/2 , when the interchanging 
inversion J12 is compensated by the permutation of the pairs of arguments. It 
gives the known Keller - Dykhne result (2). The similar fixed point in case of 
N = 3 was found by Dykhne [3j at equal phase concentrations ( Xi = x (i = 
1,2)) and the additional constraint on 03 

°3 = °o = cios, (16) 

which remains 03 unchanged under the duality transformation. The corre- 
sponding exact value for a e coincides with (2) in accordance with (15'). 

It means that at these values of partial conductivities a e is fixed and does 
not depend on xl Note that this value of a e ensures the right values in the 
limits x — > 1/2 and x — > 0, when they must be ^Jo\<J2 and 03 respectively. 

Of course, one can choose the pairs with equal concentrations by various 
ways. For example, if the phases 1 and 3 have the equal concentrations x 
and CTq = <7i<73 = cr|, then one obtains 

Ve = = (T 2 . (17) 

Analogously, when X2 = X3 and Cq = 0203 = uf, one obtains 

<7e = y/a 2 CT3 = 0"i- (18) 

The corresponding lines of equal concentrations of different phases on the con- 
centration phase diagram are shown on Fig. lb, where they coincide with the 
bisectrisses of the triangle (this triangle is a concentration space of N = 3 sys- 
tems due to the normalization condition Yli x i = 1 )■ Their unique intersection 
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point corresponds to the case of equal concentrations of all phases. It follows 
from above results that under the conditions (16)-(18) <7 e remains fixed and 
does not change along the corresponding lines! It is worthwhile to note that 
now, contrary to the N = 2 case, the exact value of <r e at the point of equal 
phase concentrations x\ = x 2 = £3 = 1/3 cannot be define from the duality 
relation in general case. It can be found only in the above fixed points and its 
value depends on the choice of fixed point. All these values are different. Due 
to the additional constraints on <jj , the form of possible exact values of a e 
admits also a value of the conductivity of the third phase. For this reason one 
can say that a structure of the fixed points and a form of corresponding exact 
values are defined by the involved interchanging inversion transformations and 
the additional constraints. 

3. General N -phase case 

We show now that the duality relation (4) for N -phase systems admits a 
generalization of these results on the N -phase systems. 

Of course, for systems with N > 3, there arc similar fixed points, when 
Xi = Xj and all other o\ (/ ^ i,j) arc equal between themselves. But this case 
is trivial, since it reduces to the 3-phase case. Fortunately, the new possibilities 
appear for N > 3. 

Let us consider firstly a case when N is even N = 2M. Then a new 
fixed point is possible, which corresponds to M pairs of phases with equal 
concentrations and equal products of the corresponding conductivities. For 
example, the fixed point 

X2i-1 = X 2 i = Ui, ol=0 2 i-\(72i (i = l,...,M) (19) 

is possible. But all <7j in the phase pairs must be different. In other words the 
conductivities of these pairs correspond to M different points on the curve (see 
Fig. la) 

2 

<7 = <7ia 2 - 

The effective conductivity has in this case the next exact value 

a e = (T = V CT 2i-io-2j (i = l, ...,M). (20) 

Note that the equal concentrations of the phase pairs y a (a = 1,...,M) can 
be arbitrary, except the normalizing condition X^fli 2y a = 1- The exact value 
(20) does not depend on y a and ensures again the correct values for a e in the 
limits y a — > 0, 1 (a = 1,...,M). Of course, the similar fixed points exist for 
other ways of the partition of N phases into the pairs with equal concentrations. 
A number of such points is equal #2M = (2M — 1)!!. The exact values of a e 
in these points have always the same general form 

u e = (t = s/oWj {i < j, i,j = 1, — , N), (21) 

where possible pairs ij correspond to the phases with equal concentrations in 
the corresponding fixed points. 
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Now we consider a case when N is odd N = 2M+1 . In this case a new fixed 
point is possible if the above conditions of the even case (19) are supplemented 
by a condition on a^u+i analogous to the Dykhne condition for N = 3 case 

°2M+1 = °0 = 02i-lCT 2 i (« = (22) 

The corresponding effective conductivity has again the exact value (20), coin- 
ciding with the exact value of system with N = 2M . Of course, the other fixed 
points related with various ways of choice of the phase pairs with equal concen- 
trations are possible. Their number is equal to #(2M+i) = (2M + 1)#2M = 
(2M + 1)!!. A general form of the exact values of a e remains the same as in 
the even case with an additional equality 

o- e = er = <r k = y/aiVj, i^j^k, i<j (i, j, k = 1, N). (23) 

where the pairs ij correspond to the phases with equal concentrations and k 
corresponds to the unpaired phases. 

It follows from (21,23) that the exact values of a e at the point, where all 
phase concentrations are equal, are not universal. 

One can show that only the fixed points of type (19,22) and their various 
permutations are admissible for all dual transformations of the form T( J, P) = 
\\JaY\Pb, where, for clarity, both products are assumed to be finite, and the 
action of the product of inversions J = JJJ a must be equal to (or compensated 
by) the product of permutations P = Y\Pb (we will not consider here possible, 
more sophisticated, forms of T). To do this it is useful to pass from partial 
conductivities a-i to their logarithms aj = lnuj, (— oo < < oo). Then the 
partial inversion transformation J„} acts on ai as a reflection a» — > A— a-i, A = 
2 In do- The inversion J ao acts on the set of variables {a} in the following way 

a t -+A-C*, (i= 1,...,A0 (24) 

The action of the product of L different inversions J = nf=i J a (°) gives anal- 
ogous result 

A L + {-l) L ai , (z = l, ...,A0 (25) 

with Al = E^i( _ l) s ^si A s = 21no-Q S '. Since all parameters of inversion 
transformations collect into one constant Al, it will be enough to consider 
only L = 1 or L = 2. In order to find any fixed point the result (25) must be 
equivalent to some permutation P of the set {a} 

{a} P = (ai p ,...,a Np ). 

Let us consider firstly a case of even L = 2. Then one obtains from (25) 

a lp =A 2 + a u (i = l,...,N) (26) 

The set of equations (26) is very restrictive. If any a, remains unchanged, than 
A2 = and all other a, also remain unchanged. It means that in this case 
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J = I. If (26) gives a permutation of any pair a; «-> d^, (i ^ k), than it follows 
from (26) that again A 2 = 0. As is known any permutation P decomposes into 
product of independent chains of simple permutations. Then, considering all 
possible chains similar to the abovementioned simplest case of pair permutation, 
one can see that equations (26) take place only when A 2 = 0. In case of odd 
L = 1 equations (25) take the form 

a ip = At-Oi, {i = l,...,N) (27) 

If (27) gives a permutation of some pair cij <-» dfc, (i 5^ k), then one obtains for 
aj,afc only one equality (since in this case both equations coincide) 

a i + a k =A 1 . (28) 

In case of any larger chains of permutation one obtains additional equalities, 
which have only trivial solution. For example, in case of 3-chain ( 1 — > 2 — > 3 — > 
1 ) they are 

o-i + a 2 = a 2 + a 3 = a x + a 3 = A±. (29) 

It follows from (29) that this is possible only if a\ — a 2 — 0,3 = a, A\ = 2a, 
but this corresponds to the one phase case. Consequently, the equations (27) 
can take place if and only if the permutation P decomposes into product of 
pair permutations (for TV = 2M ) or into product of pair permutations and one 
identical transformation (for N = 2M + 1 ) . Then one obtains in case of pair 
permutations of type a-zi-i <-> a 2 i the next equalities 

a 2i -i+a 2i = A 1 (i = l,...,M), (N = 2M), 

om-1 + a 2i = 2a 2M +i = Ai (i = l,...,M), (N = 2M + 1). (30) 

All other solutions can be obtain from (30) by various permutations. These 
solutions give exactly the fixed points (19,22) after returning to the variables 
{a}. 

4. Checking of exact results for random heterophase systems 

The existence of the fixed point hyperplanes in space of concentrations and 
a nonuniversality of the corresponding exact values at the point of all equal 
concentrations for all N > 3 can be confirmed for random heterophase systems 
in the framework of corresponding approximate schemes One needs only 

that such approximate scheme satisfies the duality relation. 

Firstly we will show how these properties look like in the EM approxima- 
tion for N -phase square lattice random resistor network (RRN) systems. This 
approximation gives the following equation for determination of a e |S] 

N 

Xi = 

. a e + a, 
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(31) 



It can be written as an algebraic equation of the N -th order 

o-f + aiaf- 1 + ... + a N ^a e + a N = 0, (32) 

where 

a-k = s k - 2s fc = -Sfc + 2s fe (fc = 1, ...,iV- 1), ajy = -s«. (33) 

Here 

N N i k _ 1 jv N 



Sk ^ 

i=l ii <i2<...<ifc_i i;=ii i=l ii<i2<---<ik »i=*i 

where all VI • , do not contain i; = i. The coefficients at are the ho- 
mogeneous functions of er, of degree k and the linear functions of concentrations 
Xi. For this reason they are related only with the corresponding symmetrical 
functions s k and their "averages" s k or s k (since there arc different ways to 
define averaged s k ). The eq. (32) must have one physical solution <r e > 0. 

It is easy to see that a e determined from the equation (32) satisfies the 
duality (and inversion) relation (4), since the equation (31) for the inverse con- 
ductivities has the next property (for all Ui ^ ) 

N 2 1 N 2 I 

E Oe — _ _ sr^ Crg/Ce - Oj _ n 

Xi ~ . — > Xi 9 . — U. 

. cr e + cr^/cr, a£/<J e + G { 

Moreover, since ajy — —sn , one can pass directly to the similar equation for / 
by dividing all its terms on sjy 

f N + a 1 f N - 1 + ... + a N ^ 1 f-l = 0, (34) 

where now the coefficients a k are expressed through projective variables Zi 

a k = a k {{z}) = a k /{s N ) k l N , k = 1, N - 1. (35) 

Firstly we consider the case of equal concentrations x% = 1 /N. Then one can 
see that the coefficients a k and / must depend only on symmetrical functions 
s k (k = l,...,N-l) 

a k = (N-2k)/N s k , f = f(h,-,s N - 1 ), (36) 

Due to the transformation rules of § k (10) the transformed / in the duality 
relation (4) has the form 

J/(si, sjv-i) = /(sjv-i, si). (37) 

Thus, one must have at the fixed point the next equalities 

s N - k = h, k = l,...,N-l. (38) 
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It is easy to check that they are satisfied at the fixed points (19), (22) and that 
the coefficients dk satisfy at these points the equations 

an-k = -a>k- (39) 

It follows from (39) that for even TV = 2M the coefficient clm — 0. Conse- 
quently, one has from (38), (39) that the polynomial in (34) can be represented 
in the factorized form 

(/-l)Pv_i(/)=0. (40) 
where a polynomial Pjv_i(/) of degree N — 1 has the next form 

M i-1 M-l M-2 M-2 

Ps-itf) = E / 2M_< E a < + / M_1 E * + ( E /')( E *). * - 2M ; 

»=i ;=o i=o ;=o 1=0 

M i M-l M-l 

^-i(/)=E/ 2M " i E a ' + (E/')(E^ JV = 2M + 1. (41) 

i=0 i=0 1=0 1=0 

Since all coefficients of Pjv-i(/) are positive, one can see from (40) and (41) 
that equations (32) and (34) have always only one physical solution at any fixed 
point (19), (22) 

f = l, a e = (sn) 1 ^ = (ijtj). (42) 

which coincides with the corresponding exact values from (21), (23). 

Now we pass to the case of pairwise equal concentrations. Unfortunately, a 
complete analysis in this case is rather complicated and, at the moment, it is 
not done in general form. For this reason we have considered only the systems 
with N = 4,5, 6. It is convenient to consider the odd and even cases separately. 
We begin with the even case N = 4, 6. Due to apparent permutation symmetry 
Sn , it will be enough to consider only fixed points (19). For N — 4,6 at the 
fixed points 

2 

Xl = x 2 = yi, X 3 = X 4 = ?/2, E = 1 A Ol °2 = 0304 = 0o, {N = 4) 

o=l 

3 

X2i-i = x 2i = yi, ^2yi = l/2, (T 2 i-i<T 2 i = ol, (i = 1,2,3), (N = 6) 
l 

one can show by direct calculation, using equalities (38), that for all yt 

a k = -a N ^ k , fe = l,2, a 2 = (N = 4), 

afe = -a N -k, k = 1, 2, 3, a 3 = (N = 6). 

Thus the polynomial in equation (34) reduces also to the factorized form 

(/-l)Pv-i(/,y) = 0, (43) 
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where the polynomials f3,s(/, y) have the form (41), but now the coefficients 
at depend on concentrations y. Due to factorization, eqs.(32),(34) have the 
required solutions of the form (20). 

In the odd case N = 2M + 1 we have checked for M = 1,2 that the 
necessary conditions = — ajy-fc, are satisfied at the fixed point (22). As a 
result the polynomial in the defining eq. (34) reduces also to the factorized form 

(/-l)iW(/,»)=0, (44) 

which means that in the odd cases the EM approximation has also the fixed point 
hyperplanes in the concentration space and that at these fixed points there are 
always the solution, which reduces at these points to the corresponding exact 
values (23). 

We are sure that the EM approximation for general N contains all fixed 
points (19,22) with the corresponding exact values. 

Now we check the obtained exact results in framework of two other approxi- 
mations applicable for some random heterophase systems. They were introduced 
firstly for two-phase systems in [5] and generalized later on N -phase systems 
in ^U]. Fortunately, they give for a e the expressions more simple than the EM 
approximation. 

We begin with the finite maximal scale average (FMSA) approximation valid 
for random systems with compact inclusions. The corresponding approximate 
expression for cr e of N -phase systems is 

N 

a e ({x},{<r}) = n<*. (45) 
i=i 

It is easy to check that (45) reproduces at the corresponding fixed points the 
exact values from (20), (21), (23) due to the normalization condition for concen- 
trations 

M 

°e = J] = °* = (sn) 1/N , (N = 2M) 

t=l 

M 
i=l 

It is interesting that in this case a e at the point of all equal concentrations 
Xi = l/N always coincide with s]^ N , (i.e. for arbitrary c^. 

Another approximate expression for a e is valid for the composite random 
heterophase models [5] . It has the following form in N -phase case ^U] 

<j e ({x},M) = V{*)(v- 1 )- 1 = f^V 7 ". (46) 

\sn-i / 

By direct substitution of the fixed points values (19), (22) into (46) one can also 
reproduce the exact values from (20), (21), (23). 
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5. Conclusion 



Using the exact duality relation, symmetry and inversion properties of 2D 
isotropic heterophase systems we have found all possible exact values of their 
effective conductivities. The obtained results take place for various heterophase 
systems (regular and nonregular as well as random), satisfying the duality re- 
lation (4), and show very interesting properties. Especially, an existence of 
hyperplanes in the concentration space with constant a e is unusual. It will be 
very desirable to check them experimentally. 
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